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When considering strangelets at finite temperature it is important to obey the constraint that any 
observed state must be color singlet. This constraint results in an increase in strangelet masses as 
calculated at fixed entropy per baryon. We use the color singlet partition function for an MIT bag, 
derived using the group theoretical projection method, to calculate strangelet masses. Mean shell 
effects are included in a liquid drop model, by using a density of states obtained from the multiple 
reflection expansion. Another important effect of the color singlet restriction, namely many orders 
of magnitude suppression of thermal nucleation of quark-gluon plasma in collisions, will also briefly 
be described. 



I. INTRODUCTION 



The (meta-) stability of strange quark matter as suggested by several authors [fL|-|4| has spurred theoretical as well as 
experimental activities aimed at answering the question of whether strange quark matter already exists in Nature (e.g, 
in the cores of neutron stars or in cosmic rays) or can be produced in a laboratory (most likely in relativistic heavy-ion 
\ collisions). It is not an easy question to answer with our present theoretical tools, since it involves non-perturbative 
• aspects of the strong interaction. Precisely for this reason most of the theoretical work has been done using different 
phenomenological models. Ultimately the question of the stability of strange quark matter will have to be answered 
experimentally/observationally, e.g, by the detection of strangelets in heavy ion collisions. Strangelet production in 
ultrarelativistic heavy-ion collisions would also constitute an unambiguous signal of quark-gluon plasma formation. 

While not making specific predictions for the outcome of an experiment here, we do try to address some of the 
points which are believed to be important for experimental strangelet searches. These are the effects of small system 
CIh' s i ze an d the fact that the produced strangelet is not born cold, but has to survive the hot environment ensuing 
the collision of two heavy nuclei. In addition to the added thermal energy, which in itself results in an increase of 
qj ' strangelet masses, there is a more subtle effect, viz. the fact that the strangelet must be color singlet. 

. We have used the MIT bag model || to describe a strangelet as an ideal Fermi gas of u, d and s quarks (and also 
antiquarks and gluons at finite temperature). We use a density of states derived using the multiple relection expansion 
H which gives corrections due to the finite system size. This is the liquid drop model of strangelets, similar to the 
liquid drop model of nuclei. The overall effects of the shells are well described by the liquid drop model, whereas 
the detailed effects of the shells, such as shell closures, which are likely to be important at very low baryon number 
- ■ are not revealed in the liquid drop model. However, recent results indicate that shell effects are washed out at 
elevated temperatures, making shell model calculations and liquid drop calculations effectively equivalent. The main 
advantage of the liquid drop model over the shell model is that many quantities can be evaluated analytically. 
Some of the results presented here are described in Ref. H, where only massless quarks were treated. 



II. STRANGELETS AT FINITE TEMPERATURE 



We first consider strangelets at non-zero temperatures without inclusion of the color-singlet constraint. In this case 
we can obtain all relevant information from the grand canonical partition function Z of a noninteracting gas of quarks, 
antiquarks, and gluons at temperature T — fi^ 1 in an MIT bag of volume V 

In Z = ln Z i ~ PBV. (1) 

i 

Here Zi is the partition function for particle species i and B is the bag constant, which is the energy density of the 
perturbative vacuum inside the bag. Quarks, antiquarks, and gluons contribute terms of the form 
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lnZi = ±gi / dk p(k) In {1 ± exp [- (e(k) - fii)}}, 



(2) 



where the upper sign is for fermions, and the lower for bosons, gi is the statistical weight due to the spin (helicity) 
and color degrees of freedom, e(k) is the energy of a particle with momentum k, and fi is the chemical potential. We 
use units in which h = kg = c = 1. The density of states p(k) is a smooth function of k calculated in the framework 
of the multiple reflection expansion o . It has the general form 
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where S is the surface area of the bag, and C = J s dS{\/R\ + I/R2) is the extrinsic curvature of the bag surface (Ri 
and i?2 are the principal radii of curvature) . The functions fs and fc depend on the equations of motion and on the 
boundary conditions, but not on the geometry of the confining surface, which is only assumed to be smooth. The 
surface coefficient for a quark field in a cavity with confining MIT bag boundary conditions is [s|JlO|] 
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Note that this expression vanishes in the limit m — > 0, so that massless quarks do not contribute to the surface 
tension. So for massless (u and d) quarks the lowest order correction to the density of states is the curvature term. 
The curvature coefficient fc has not been calculated using the multiple reflection expansion in the general case of 
massive quarks, but for massless quarks with MIT bag boundary conditions the result is fc = — 1/247T 2 . In the 
limit of infinite mass the quark can be considered non-relativistic and the field equation reduces to the wave equation. 
In this case the boundary condition of the MIT bag is the well known Dirichlet boundary condition studied by Balian 
and Bloch B who derive the result fc — 1/12-7T 2 . It was shown by Madsen pf that the expression 
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which has the above special cases as limits, is in very good agreement with shell model calculations (see also Fig. [|). 
We will therefore adop t (H ) as the curvature coefficient for quarks. For gluons the relevant expressions were given by 
Mardor and Svetitsky flflf 



fs = 



J( ' 6tt 2 ' 



(6) 



Just as massless quarks, gluons do not contribute to the surface tension. 

It is possible to calculate the thermodynamical potential £1 = —T In Z without resorting to numerical evaluation in 
the two limits: T — » and m — * 0. Since we are mainly interested in the finite temperature regime we will give the 
result for J\f q massless quark flavors, including antiquarks and gluons, and assuming a common chemical potential \x q 
of all quarks (antiquarks have fig = —fi q ) 



Q(T,p q ,V,C) 
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In the following we consider only spherical systems, where V ~ 47ri? 3 /3 and C = 8ttR are given in terms of the radius 
R. When fixing the geometry in this manner there is only one independent "shape variable" which we will take to 
be the volume. A strangelet is in mechanical equilibrium when <9f2/ 'dV\ ll t = 0. The chemical potential and the 
temperature are determined by requiring the baryon number A and the entropy per baryon S/A to be fixed. So to 
find the equilibrium values of T, p q and V we simultaneously solve the equations 



dV 



^77 =0, 



on 

dflg / rp y 



3A, 



1 (dfl 
A (df 



The energy (mass) E of the strangelet can then be calculated by use of the identity 
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(8) 
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E = il + TS + 3/iqA. 



(9) 



As long as only massless quarks are considered this is equivalent to the expression E = ABV in equilibrium. This is 
because massless particles satisfy the relativistic equation of state E — 3PV , where P is the pressure. In equilibrium 
the dynamic pressure of the particles is exactly balanced by the bag constant B = P, so that adding the vacuum 
energy BV gives E = 3BV + BV = ABV. Including a massive strange quark this no longer holds. Plots of the energy 
per baryon as a function of baryon number obtained by solving Eqs. (g) are shown as full lines in Fig. [p. 
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FIG. 1. The energy per baryon number as a function of baryon number. Full curves show results without the color singlet 
constraint — dashed curves with it. Left: curves for fixed entropy per baryon (from top to bottom: S/A = 10,5,1,0) with three 
massless flavors. Right: Equivalent results for a strange quark mass of 150 MeV. The two sets of curves correspond to S/A = 
and 5. The jagged curve is a shell model calculation. In both panels B 1 ^ 4 = 145 MeV. 



III. COLOR SINGLET STRANGELETS 

The partition function given by Eqs. ([!]) and (||) is in fact the grand canonical partition function for a colorless 
assembly of quarks, antiquarks, and gluons in an MIT bag. In order to force a non singlet color configuration it 
is necessary to introduce chemical potentials related to the color degree of freedom of the quarks and gluons. Two 
chemical potentials, corresponding to the two conserved charges (Cartan operators) of SU(3), are needed. Now we 
need to put in a sum over the different color states of quarks and gluons in the partition function, thus replacing 
Eq. (g) with 

lnZ i =±g' i y2 i dkp{k)\n\l±e- ^-^-^ c ^-^ c ^\. (10) 
c Jo 1 J 

The statistical weight g\ no longer includes the color degree of freedom, which is now explicitly summed over. C% and 
C 2 are the eigenvalues of the Cartan operators, which for quarks are in the fundamental representation and for gluons 
in the adjoint representation. Using the group theoretical projection method developed by Redlich and Turko [[12] fl4| , 
it is possible to construct a canonical partition function (with respect to color) from the grand canonical ^(71,72) 
constructed using the prescription just outlined. The resulting color-singlet canonical partition function is 

Zc.s. = [ dd 1 d9 2 M(9 1 ,6 2 )Z{iT9 1 ,iT0 2 ), (11) 

J — IT 

where M(9i,6 2 )d9id9 2 is the Haar measure of SU(3). Note that the chemical potentials ji — iTOi are now purely 
imaginary. The baryon number is still treated in a grand canonical ensemble, so Zc.s. is a mixed canonical/grand 
canonical ensemble. 
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In the thermodynamic limit (V — > oo) the two ensembles are equivalent, i.e., they give identical mean values — 
differing only with respect to fluctuations. For a finite system this is no longer the case, and since color neutrality is 
exact at all times the canonical prescription is the proper one. 

For massless quarks it is possible to evaluate Zq.s. analytically in a saddle-point approximation, valid at high T 
and/or [i q . This was done by Elze and Greiner |lq] , who for Af q massless quark flavors obtained the result 



Z c .s. (T, V, n g ) = n(T, V, » q )Z(T, V, ft q ), 



(12) 



where II is given by 
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Note that II vanishes in the limit V — ► oo. The limit T — ► cannot be taken in the above expression for II, since the 
saddle-point approximation breaks down for low T. But for T — > the effect of exact color singletness also disappears 
since at T = it is always possible to construct a color singlet state for 3A quarks occupying the lowest energy levels 
available. For massive quarks we evaluate the partition function numerically. 

The effect of the color singlet constraint on the energy per baryon is seen in Fig. [l], where dashed lines include the 
color singlet constraint. There is a sizeable effect regardless of whether one includes a massive strange quark or not. 
The main effect of the strange quark mass is to shift curves towards higher E/A. Curves for zero strange quark mass 
are calculated using the saddle-point approximation, which has been shown to be a good approximation in this case 
(see Ref. §). 



IV. COLOR SINGLET SUPPRESSED NUCLEATION 



As another example of the effect of color singletness we consider the nucleation of quark matter droplets from 
hadronic matter using standard homogenous nucleation theory. The low baryon number density region of the phase 
diagram is where the transition to a quark-gluon plasma would happen in RHIC and LHC experiments. In the region 
of the phase diagram characterized by high baryon number density and moderate temperature this could be relevant 
for the nucleation of quark matter droplets in AGS experiments and in the core of neutron stars. In the case of neutron 
stars the time scales are probably sufficently large for homogenous nucleation theory to be applicable, whereas this is 
questionable in the case of quark-gluon plasma formation in heavy ion collisions. An estimate for the nucleation rate 
is 

Tl «T 4 cxp(-A J F/T), (14) 

where AF is the height of the free energy barrier that has to be surmounted by the nucleated bubble. For the 
prcfactor we have chosen the dimensional estimate T 4 . One could consider improvements to this estimate but since 
the relative suppression due to color singlet effects are dominated by the exponential, this is unlikely to change the 
results significantly. 

The free energy barrier is given by the pressure difference between the hadronic phase and the plasma phase (times 
volume) plus additional contributions from surface tension and the curvature term. 

AF = (P h - P q )V + <tS + jC. (15) 

Note that any contribution from the hadronic phase, i.e., either a non-zero pressure or a surface tension, will increase 
the barrier height. We will therefore ignore the hadron gas contribution, since it will only augment the effect due to 
color singletness (for a more elaborate discussion of this point, see Ref. P6[). With this simplification the free energy 
barrier is simply AF = Q, where Q is given by Eq. (g) with N q = 2. The color singlet constraint is taken into account 
by using £1 = — TlnZc.s. as given by Eqs. (|l2| ) and (|l3|). The barrier is shown for a particular choice of parameters 
in Fig. |. 
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FIG. 2. Left: Free energy barrier for T = 170 MeV and [s, q = 0. The upper curve is with the singlet constraint, the lower is 
not. Right: Nucleation rate for fi q = (solid), 100 (dotted), 300 (dashed), 400 (dot-dashed), and 500 MeV (dash-triple dot). 
Lower curves with singlet constraint, upper without. In both panels B 1 ^ 4 = 200 MeV. 

Inserting the value of AF for a critical bubble in the expression ( |l4| ) for a range in temperature and chemical 
potential gives a plot like the one shown in the right panel of Fig. 0, where an extra constraint, namely requiring the 
droplet to have zero momentum (lcfl , is also included. The effect of the color singlet constraint, which is the dominant 
of the two constraints, is seen to be quite dramatic. It causes a suppression of the plasma formation rate by four to 
five orders of magnitude. 



V. CONCLUSIONS 



We have shown that the effect of exact color singletness plays an important role in strangelets — increasing masses 
as calculated at fixed entropy per baryon — as well as in the phase transition between hadronic and quark matter, 
where it causes a suppression of quark gluon plasma nucleation. 
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